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We explore molecular states of two open heavy-quark mesons {Qq)-(qQ) in a quark-based model 
in terms of a four-body non-relativistic Hamiltonian with pairwise effective interactions. Molecular 
states are found in the combinations of {D,D*,B,B*} with {D, D* , B, B*}, including a weakly- 
bound DD* state near the threshold which may be qualitatively identified as the 3872 state observed 
recently by the Belle Collaboration. 
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' The narrow 3872 MeV state recently discovered by the Belle Collaboration [1] , and subsequently confirmed by the 
, CDFII Collaboration [2], has generated a great deal of interest. The state contains a cc pair and was observed in 
' the exclusive decay of —> K^tt^tt^ J/tp. Possible explanations of the state include a conventional charmonium 
state, hybrid state, and multi-quark molecular state [3-12]. Direct experimental search using e'^e" annihilation by 
the BES Collaboration in the BEPC indicates that this state is unlikely to be a vector 1 state [13]. The 3872 state 
J> ] lies about 60 MeV higher than the weighted average of the conventional ^'(IZ?) charmonium state predicted by the 
OO • Cornell potential [14] and the Buchmiiller-Tye potential [15]. The identification of the 3872 state as ^(^I?2) would 
imply too large a splitting between ^'(■^Z?2) and '^(^Di), as the latter state has been identified with the 3770 state. 
^— ^ Furthermore, the 3872 state is seen in the tt+tt" transition to J/ip and only an upper bound on El transition to 
I ' the xi'^Pj) state can be determined experimentally, whereas the ^(^1)2) state would have yielded a large partial El 
I transition width, about 5 times greater than the width for the tt+tt^ transition [16]. On the other hand, the 3872 
' state lies very near the D*^D'^ threshold. The proximity of this state near the threshold and the difficulties with the 
^(■^^2) description led the Belle Collaboration to suggest that the observed state may be a multi-quark 'molecular 
state' studied in previous theoretical investigations [3-11, 17-38]. 
r^" We shall be interested in molecular states of the type {Qq)-{qQ) where Q is a charm or bottom quark and g is a 
Oh! u 01 d quark. We shall often use the simplified nomenclature of "M and M*" or MM* to denote mesons M and 
^ ' M* or their charge conjugate M* and M. Molecular states formed by D and D* were predicted many years ago by 
Tornqvist [23] who suggested that open heavy-quark pairs form deuteron-like meson- meson states, 'deusons', because 
of the strong pion exchange interaction. In the heavy meson sector, Tornqvist estimated that one-pion exchange 
alone is strong enough to form deuteron-like composites of BB* and B*B* bound by about 50 MeV, and composites 
of DD* and D*D* bound by pion exchange alone are expected near the threshold. Manohar and Weise have also 
d ' studied another type of multi-quark systems, the (QQqq) hadronic states, in the limit where the quark mass mg of Q 
goes to infinity. They noted that there is a long-range binding due to the one-pion exchange between ground state Qq 
mesons. They suggested that for two open bottom mesons, this long-range interaction may be sufficient to produce 
a weakly-bound two- meson state [24] . Ericson and Karl discussed the physics of pion exchange and the role of tensor 
forces in forming such hadronic molecules [25]. 

The pion-exchange model presents a reasonable description of the long-range attraction between D and D* . How- 
ever, as the state energy depends on the strength of the potential at short distances, the results of the pion-exchange 
model depends on how the long-range potential is regularized at short distances. The results may also be affected 
when exchanges of more mesons are included. 

To describe the short-distance behavior, it is desirable to study multi-quark molecular states in a quark-based model. 
The phenomenological interaction between a quark and an antiquark in a single meson is reasonably well known 
[14, 15, 39-44]. However, there are considerable differences and uncertainties in the description of the interaction 
between constituents in a baryon, as it has been given in terms of one-gluon exchange (OGE) interactions involving 
color operators of the type A(i) • A(j) [26-38], or alternatively in terms of a one-boson meson exchange (OBE) 
interactions involving fiavor operators of the type • [45-49]. Here, X{i) is the generator of the SU(3)coior 
group for particle i and T(i) is the generator of the SU(2)fiavor group. The interquark interactions in the OGE and 
OBE models have very different color structures. 

For the interaction between constituents of different mesons, Lipkin and Greenberg pointed out pathological prob- 
lems involving the use of confining one-gluon exchange interactions in a quark-based model because they do not respect 
local color gauge invariance [50, 51]. The underlying difficulty arises because the quark-based one-gluon exchange 
model represents a truncation of the basis states with the neglect of the dynamics of explicit gluon degrees of freedom. 
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In spite of these defects, the quark-based one-gluon exchange model with the interchange of quarks has been used 
successfully to study the interaction of mesons at short distances [26, 40-44] when the dynamics of the gluon degrees 

of freedom can bo neglected. 

In a recent quark-based model study of molecular states using the one-gluon exchange model with constituent 
interchanges, Swanson found that the one-gluon exchange interaction couples DD* with u}{J/tp), but the interaction 
is not attractive enough to lead to a boimd molecular state [11] . Swanson then studied molecular states with the meson- 
based onc-pion exchange model (as in Tornqvist [23] ) and added the quark-based one-gluon exchange interaction with 
constituent-interchanges. The probability of the mixing of ui{J/'ip) with DD* depends on the one-pion exchange 
potential cut-off parameter A, ranging from zero for the cut-off value of A used by Tornqvist without the one-gluon 
exchange interaction [23] to a maximum mixing of 17% at large values of A. In such a treatment in which the 
interaction between quarks is superimposed on the onc-pion exchange interaction between mesons, the position of 
the molecular state relative to the DD* state and the degree of ui{J/tp) mixing depends on how the long-range pion 
exchange potential is regularized at short distances and the result may also be affected when exchanges of more mesons 
arc included. 

We search for a description of the heavy-quark molecular states within a completely quark-based model. We 
note that the physical conditions appropriate for the molecular state of our interest have important effects on the 
interaction between the constituents in different mesons. The molecular state at the center of our attention has a 
binding energy of a few MeV. For such a weakly-bound molecular state, the average separation R between the heavy 
mesons is considerably greater than the average radius of heavy quark mesons a (See Section III for more detail), as 
was already recognized by Close and Page, Voloshin, and Braaten and Kusunoki [4, 6, 10]. At these large distances 
with R >> a, the probability for the recombination of a quark of one meson with the antiquark of the other meson 
after a gluon exchange is highly suppressed as their quark- antiquark separation R is much greater than their average 
natural meson radius a. When R » a, far more likely after the exchange of a single gluon is the occurrence of the 
exchange of an additional gluon between quarks. This two-gluon emission leads to the color van der Waals interaction 
proportional to the inverse-power of R, as shown by Applequist et al., Pcskin et al., Bhanot et al., Lipkin, Grccnberg, 
and many other authors [50-55, 60]. The color van der Waals interaction can be equivalently represented in terms 
of effective charges of quarks and antiquarks in a QED-type interaction. There is however no experimental evidence 
for the color van der Waals interaction of the inverse-power type between separated hadrons. As pointed out by 
Lipkin and Greenberg [50, 51], it is necessary to include additional dynamics of gluons in such a quark-based model. 
Accordingly, due to the breaking of the gluonic string at large distance, the color van der Waals interaction should 
be screened by introducing a screening mass /x, which will modify the inverse-power color van der Waals interaction 
at large distances. Appropriate for the peculiar physical conditions of weakly bound molecular states, the concepts of 
effective charges and the screening mass arc incorporated into the effective interquark potential which is subsequently 
used to explore meson-meson molecular states. 

Nucleus-nucleus molecular states have been observed previously in the collision of light nuclei near the Coulomb 
barrier [56 58]. Effective nucleon-nucleon interactions have been successfully applied to obtain the potential between 
two nuclei to study their reactions [58]. In a similar way, effective interquark potentials can be used to obtain the 
potential between mesons. Just as with the nucleus-nucleus potential, the meson-meson potential can be evaluated 
as a sum of a direct potential and a polarization potential. The direct potential arises from the interaction of a 
constituent of one meson with a constituent of another meson, and the polarization potential arises when one meson 
polarizes the other meson in its vicinity. Using Gaussian wave functions for the mesons, the direct potential and the 
polarization potential can be obtained analytically. The knowledge of the potential between the mesons then allows 
one to determine the eigenstate of the four-quark system. 

Using the; c;ffectivc! interaction, we find many weakly-bound moleciilar states in {Qq)-{qQ) systems with charm and 
bottom quarks, including a weakly-bound D~^D*~ state near the threshold which may be qualitatively identified as 
the 3872 state. For the bottom meson {Qq)-{qQ) pairs, there are weakly-bound 28 states near the threshold, in 
addition to IS states with binding energies of about 150 MeV. 

Our exploratory study using effective charges and effective interactions in a non-relativistic quark model differs from 
those of the pion-exchange models and the multi-quark models with explicit single- gluon or single-boson exchanges. 
The results of molecular states obtained from different models can naturally be different. A careful comparison of 
these different results with experimental data will be useful to determine the importance of various mechanisms that 
are present in generating the molecular states. 

This paper is organized as follows. In Section II, we start with the four-quark Hamiltonian and partition the 
Hamiltonian into an unperturbed part of two mesons and a residue interaction. The eigenvalue equation for the 
four-quark system is reduced to a Sclirodingcr equation for the relative motion of the two mesons in a meson-meson 
potential. In Section III, we present justifications for the introduction of the effective charges and the screening mass 
in the effective interquark interaction. In Section IV, we present our method for the evaluation of the molecular 
potential. Section V presents the results of the meson-meson potential and eigenvalues for four-quark states. We find 
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molecular states in the combinations of {D,D*,B,B*} with {D,D* ,B,B*}. In Section VI, we discuss the origin of 
the short-distance attraction in {Qq)-{qQ) systems and contrast them with the short-distance repulsion in {Qq)-{Qq) 
systems. In Section VII we present our conclusions and discussions. 

II. HAMILTONIAN FOR THE FOUR-QUARK SYSTEM 

To study the {Qq)-{qQ) system, we label constituents Q, q, q, and Q as particles 1, 2, 3, and 4 respectively and 
describe the four-quark system with a non-relativistic Hamiltonian 

j=l j=l k>j j=l 

in which particle j has a momentum pj and a rest mass rrij. The pairwise interaction Vjk{rjk) between particle j and 
particle k depends on the relative coordinate between them, 

Tjk = Vj - Tk- (2) 

We introduce the two-body momentum 

Pjk=Pj+Pk, (3) 

and the two-body internal relative momentum 

Pjk = fkPj - fjPk' (4) 

where 

fk = mk/nijk, (5) 

mjk = mj+mk. (6) 

There are many ways to partition the total Hamiltonian H into an unperturbed two-meson part and a residue 

interaction Vj . Wc partition the Hamiltonian so that in the lowest order the state in question can be described by a 
state of the unperturbed Hamiltonian. We can, for example, choose to partition H into an unperturbed Hamiltonian 
hi2 + of mesons A{12) and B(34) according to 

H = ^ + ^ + Vi + hi2 + hs4, (7) 
2mi2 2m34 

Vi = Vii{ru) + Fi3(ri3) + F23(r23) + V'24(r24), (8) 

hjk = ^+Vjk{rj-rk) + mjk for (jfc) = A(12) and 5(34), (9) 

Z/djk 

where jjjk = rrijmk/rnjk- The eigenvalues of the Hamiltonians hi2 and ^34 can be solved separately to obtain the 
bound state wave functions and masses Mjk{y) of mesons A and B, 

hjk\{jk)u) = + mjk\\{jk),) = Mjk{iy)\ijk),). (10) 

The four-body Hamiltonian becomes 

H=P^ + p^ + Vi + M,2{u)+Mu{iy). (11) 

2mi2 2m34 

The above non-relativistic approximation is obtained from relativistic results by neglecting terms of order Cjk/Mjk 
and higher [59] . In order to satisiy the boundary condition at large separations for which Vj approaches zero, we need 
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to include some of these higher-order terms and modify TO12 and 77134 in the above equation to Mi2{v) and M^i{v') 
so that the Hamiltonian becomes 



H = 



+ 



34 



2Mi2(z^) 2M34(!/') 



- + y/ + Mi2(!^)+M34(l/')- 



(12) 



The above Hamiltonian then describes properly the asymptotic behavior at large separations of the two mesons. 

To solve for the cigcnstatcs, wc use the two-body meson states |A^(12)i?^' (34)) as basis states. We consider a 
four-quark system in which the mesons would be in the \A{12)^B{2>^)^i) state if there were no residue interactions 
Vi. When we include the residue interactions V/ as a perturbation, the eigenfunction of H becomes 



where 

r = R12 — R34, 
and Rjk is the center-of-mass coordinate of rrij and m^, 

Rjk = fjrj + fkVk- 

The spatial coordinates of the four-particle system and r, rjk, and Rjk are shown in Fig. 1. 
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FIG. 1: The spatial coordinates of the four-particle system. 

The wave function ^/'(r) describes the relative motion between the two mesons and the symbol Y^'xy to indicate 
that the sum is over all meson states except \Ai,Byi). The eigenvalue equation is 



(16) 



Working in the center-of-mass frame and taking the scalar product of the above equation with \AvBi,'), we obtain 
the Schrodinger equation for relative motion of mesons A^{12) and i?y'(34), 



2/x 



AB 



+ V{r) \ ijj{r) = eip{r), 



where p is the relative momentum 



P = 



and hab is the reduced mass of the two mesons 



IJ-AB = 



M^2{v)M^i{v') 



(17) 



(18) 



(19) 
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The meson- meson potential V{r) is given by 



We shall call the first term the direct potential 

VdUr) - (A,(12)B,.(34)|V/(r,ri2,r34)|^.(12)S,y(34)). (21) 

It arises by the direct interaction of the quark matter densities of the projectile meson with the quark matter density 
of the target meson mediated by the residue interaction Vj. We shall call the second term the polarization potential 

It is always attractive and it arises from the excitation of the colliding mesons into intermediate states due to the 
residue interaction Vj, as in the polarization of a meson in the vicinity of another meson. When the potential V{r) is 
determined and the eigenvalue e of the Schrodinger equation [Eq. (17)] is obtained, the mass of the four-quark system 
will be given by M = Mi2(z^) + M34(i^') + e. 



III. EFFECTIVE INTERACTIONS 



We shall consider the lowest energy state for which the unperturbed mesons are in their ground states with u = 
u' = 0. The energy of the four-quark system depends on the interquark interaction. The low-energy properties of 
isolated mesons can be described reasonably well by a quark and an antiquark interacting with a phenomenological 
confining onc-gluon exchange interaction [14, 15, 39 44]. When generalized to the multi-quark interactions involving 
two quark-antiquark pairs, the one-gluon exchange interaction contains a color operator of the type F(i) ■ F[j), 
where F{i) is A(i)/2 for a quark and is (— A^(i)/2) for an antiquark. It was realized by Lipkin and Greenberg that 
such a quark-based one-gluon exchange model violates local color gauge invariancc because of the neglect of the 
gluon dynamical degrees of freedom [50, 51]. Nevertheless, the quark-based one-gluon exchange model has been used 
successfully to study the interaction between mesons at short distances [26, 40-44] when the gluon dynamical degrees 
of freedom can be neglected. 
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FIG. 2: Two different ways to neutralize the color after exchanging a gluon at time to' (a) Color can be neutralized by the 
C interacting with the C to form a J/tp, and the q interacting with the q to form a w or p . (6) Color can be alternatively 
neutralized by exchanging another gluon at ti. The multi-gluon exchanges in the initial and final states represent schematically 
non-perturbative QCD interactions leading to bound meson states. R is the separation between mesons and a is the average 
meson radius. 



When one applies the one-gluon exchange model to study the interaction between mesons, the two mesons become 
color octet states after a gluon is exchanged between them. The colors of the two mesons can be neutralized when a 
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quark of one meson interacts with the antiquark of the other meson to form a bound color-singlet state (Fig. 2a), as 
described by the quark interchange model of Barnes and Swanson [40, 41]. When this quark-based model is applied 
to study heavy-quark molecular states, Swanson found that the onc-ghion exchange interaction couples DD* with 
u){J/tp), but the interaction is not attractive enough to lead to a bound molecular state [11]. Swanson then studied 
molecular states with the meson-based one-pion exchange model and added the quark-based one-gluon exchange 
interaction with constituent-interchanges. He found that the probability of the mixing of oj{J /ip) with DD* depends 
on the onc-pion exchange potential cut-off parameter A, ranging from zero for the cut-off value of A used by Tornqvist 
without the one-gluon exchange interaction [23] to a maximum mixing of 17% at large values of A. 

To study heavy-quark molecular states within a completely quark-based model, we note that the physical conditions 
associated with the molecular state of our interest have important effects on the interquark interaction. We consider, 
for example, the molecular state of D^D*^ to be a possible description of the 3872 MeV state. It is then a weakly- 
bound molecular state with a binding energy e given by e = [3.872GeV - M{D+) - M{D*-)] = 7.4 MeV. The 
average separation R between £)+ and D*~ is of the order of i? ~ h/ \/2rrh^ where mred is the reduced mass, 
M{D+)M {D*~)/{M {D+) + M {D*~)). For the binding energy of 7.4 MeV, we can estimate the average separation R 
between D^ and D*~ to be about 1.6 fm. This separation R is much greater than the average radius a of D^ and D*~ , 
whose root-mean-square radius is about 0.3 £m as given by half of the calculated root-mean-square quark- antiquark 
separation listed in the Appendix of [42] . Our molecular state is therefore characterized by an average meson-meson 
separation R much greater than the average radius a of the mesons. This peculiar feature of i? >> a has already 
been recognized by Close and Page, Voloshin, and Braaten and Kusunoki [4, 6, 10]. 

After a gluon is exchange between £)+ and D*~ at Iq, the C and C are separated by a distance of the order of 
the meson-meson separation R, as depicted in Fig. 2a. On the other hand, a bound C and C will have a radius of 
order a. When R >> a, there is a mis- match between R and a. Because of this mis-match, the probability for the 
C to interact with the C to form a bound color-singlet CC after the exchange of a single gluon, as depicted in Fig. 
2a, is highly suppressed. The estimate of Swanson[ll] on the mixing of the quark-interchange component of Fig. 2a 
(using a superposition of quark-based and meson-based models) gives a mixing probability ranging from zero percent 
to a maximum of 17 percent. For R >> a, this one-gluon exchange contribution as represented by Fig 2a is not the 
dominant process. Within a quark-based model, far more likely after the exchange of a single gluon is the exchange 
of an additional gluon from one meson to the other meson to neutralize the octet colors, as depicted in Fig. 2b. For 
our case of i? >> a, it is reasonable to consider only Fig. 2b as the dominant contribution in the present exploratory 
work. The contribution of Fig. 2a can be included in future refinements. 

The time scale for a meson to be in a color-octet state after emitting a single gluon, which can be called the color- 
octet persistence time, is of order h/{as/a) [53]. The time scale for molecular motion is of order H/e. Therefore, the 
color-octet persistence time is much shortcir than the time for molecular motion. As a conseqiience, mesons undergoing 
molecular motion are predominantly in their color-singlet state. We can invoke arguments analogous to those used 
previously by Applequist et al. [52] and Peskin and Bhanot et al. [53-55] to treat the many gluons emitted from a 
meson to the other meson during the short color-octet persistence period as grouping together into clusters, and each 
multi-gluon cluster summing to a total color-singlet. 

Peskin evaluated the set of 10 two-gluon emission diagrams which occur in the short color-octet persistence period 
[53]. He showed that in order to obtain gauge invariant results, the gluons originating from a meson are emitted not 
only from the quark and the antiquark of the meson, but also from the gluons exchanged between the quark and 
the antiquark within the meson. After the QCD two-gluon emission diagrams from a meson have been summed, the 
coupling of gluons to a static meson source is similar to the familiar coupling of QED photons through simple QED 
electric dipole interactions. The interaction between mesons can be approximated by a dipole-dipole interaction as in 
the van der Waals intc;raction between electric dipoles [53 55] . The only new feature is that the quark and antiquark 
of the meson is bound into a dipole pair with effective charges for the quark and the antiquark. These similarities will 
allow us to introduce the effective charges in a QED-type interaction to describe the interaction between constituents 
in different mesons. 

We can examine the results of Peskin and Bhanot to motivate our introduction of the effective charges. In QCD, 
after the 10 two-gluon emission diagrams have been added together, Peskin [53] and Bhanot and his collaborators 
[54, 55] obtained the gauge-invariant result that the square of the color-electric field strength at a distance r changes 
from the r"'' behavior at very large distances to a r^® behavior at intermediate distance, in the same way as the 
Casimir-Polder effect in atomic physics [61]. At a distance r with a/ag > r > a, in the vicinity of a color-singlet 
quark- antiquark meson represented by a Wilson loop, the square of the color-electric field strength is [53-55] 

(^^W) = ^(l + 3cos^^), (23) 

where C2 = (iV^ — 1)/2N is the Casimir eigenvalue for the fundamental representation of the SU(N)coior group, 
as = g'^ l^'K, g is the strong interaction coupling constant, and 6 is the angle between a and r. On the other hand. 
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the square of the electric field strength in QED at a distance r, with a/a > r > a, in the vicinity of an electric charge 
dipole represented by a Wilson loop, is [55] (see also [62]) 



aa 



2 



(^^QED(r-)) = ^(l+3cos^e), (24) 

where a = e'^/in is the fine-structure constant. Therefore, if the quark has an effective charge Cq = \fC2 (in unit 
of g) and the antiquark an equal and opposite effective charge Cq = —\fC^ (in unit of g), and each effective charge 
generates a color electric field strength as as in QED, then the square of the total color electric field strength in 
the vicinity of a meson will be the same as given in Eq. (23) obtained by summing up two-gluon emission diagrams in 
QCD. Because of the similarities between Eqs. (23) and (24), effective charges in a QED-type interaction is a simple 
way to represent multi-quark QCD interactions in the vicinity of a meson. At short distances, the effective charges in 
a QED-type interaction also leads to the proper QCD color-Coulomb interaction between a quark and an antiquark 
in a color-singlet state 

VM = (25) 

The interaction between a quark and an antiquark contains not only of the color-Coulomb component but also 
the linear confining potential. We need to modify the color-Coulomb interaction to include both the color-Coulomb 
potential and the linear confining potential. Recent lattice gauge calculations show that the nonperturbative in- 
terquark potential between a quark and an antiquark in different representations is proportional to the eigenvalue 

of the quadratic Casimir operator C2 [63]. Thus, there is the Casimir scaling in the interaction of a quark and an 
antiquark and the concept of the effective charge is also applicable when the potential is generalized to include both 
the color-Coulomb and the confining linear components. 

Color interaction given by Eq. (23) is the well-known van der Waals interactions arising from repeated applications 
of the one-gluon exchange interaction for extended and separated mesons [50-55, 60]. There is however no experi- 
mental evidence for the inverse-power color van der Waals interaction between separated mesons. This experimental 
contradiction led Lipkin and Greenberg [50] to note that the quark degree of freedom with only one gluon-exchange 
interaction for a multi-quark system is incomplete. The gluon in non-Abelian QCD with its color charges and nonlin- 
ear self-interaction is very different from the neutral photon. The gluonic degree of freedom is needed in multi-quark 
interactions [50]. One can represent the gluonic degree of freedom in terms of a string, and the important effect 
of the dynamical string is the breaking of the string at large distances, resulting in the screening of the long-range 
interaction. Accordingly, we should modify the phenomenological "color- Coulomb plus linear" interaction further by 
introducing a screening mass /x, because we intend to use the effective interaction also for distances larger than 1 fm 
for our molecular state. 

When Eq. (25) is written as Vcif) = ^C2Vc{r) by separating out the color factor {—C2), the basic color-Coulomb 
interaction is vdr) = Og/r and the corresponding basic linear confining interaction is vunir) — —Zhr/A. They can 
be represented in the momentum space by Vc{k) = Awas/k^ and ■uiin(fc) = Gnb/k'^. The effect of screening can be 
introduced by replacing with k^ + M^, leading to an interaction in momentum space given by [64] 



^J(k) = 



k2 + ^ (k2 + ^2)2 



(26) 



From the interaction in momentum space, we can obtain the interaction in the the configuration space. The screened 
color-Coulomb interaction becomes the Yukawa interaction, and the screened linear potential becomes the exponential 
interaction. One also needs to include the hyperfinc spin-spin interaction. As the spin-spin interaction is short-ranged, 
the effect of screening on the spin-spin interaction can be neglected. 

The effective interaction between quark and antiquark particles j and k in the four-body system is then given by 



where 



and 



Vjk{r jk) = CjCkv{r jk) (27) 
\/C2 if j is a quark. 



—\fC2 if j is an antiquark. 



v{rjk) = + -r^ ' - o ^3 ■ e " "^iK 29 
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The central potential reduces to the usual Cornell- type potential in the limit when /j, 0. With this effective 
interaction, a quark and an antiquark in the four-body molecular state {Qq)-{qQ) system are subject to the same 

phenomcnological interaction and can lead to the proper meson bound states. 
The residue interaction between A{12) and .6(34) is then 

2 4 

Vr = ^ ^ CjCkv{rjk)- (30) 

j=l fe=3 

By making use of earlier studies of potential parameters and a running coupling constant [42], we use an effective 
interaction with the following set of parameters for our present work 

127r 

"-(^^^ = (33^2n,)ln(^ + QVi?-) ' ^ = ^ = '-'^''^^^ 

6 = 0.335 GeV^ ^ = 0.28 GeV, ct = 0.897 GeV, (31) 
rriu = md = 0.334 GeV, = 1.87 GeV, rnt = 5.18 GeV. 

The string tension coefficient h = 0.335 GeV^ and the screening mass ji = 0.28 GeV were previously found to give 
reasonable descriptions of the charmonium masses [65]. The value of the screening mass gives a screening length of 
1/ jjL ~ 0.7 fm which is consistent with the minimum string-breaking distance Lmin, (Eq. (5.11) of [66]), based on the 
Schwinger mechanism [67] of pair production, 

Lmin = 2m„,d/K = 0.67fm, (32) 

where we have used a string tension coefficient k of 1 GeV/fm. The large value of b in Eq. (31) arises because the 
effective string tension coefficient 6e~'"'''= is not constant and it varies with rjk- In Eq. (31), we identify Q as the sum 
of the rest masses of the interacting quarks. The parameterization of the phenomcnological strong coupling constant 
as and the spin-spin interaction parameter a come from earlier calculations of meson masses in [42]. The set of 
parameters in Eq. (31) leads to a good description of the masses of D, D*, J/'ip, B, B* , and T mesons. 



IV. EVALUATION OF THE MESON-MESON POTENTIAL V{r) 

We represent the wave function for the coordinate of the quark Q relative to the antiquark q in the meson A{12) 
by a Gaussian wave function of the form 



"^^-^"-'^^ = V (^TTT)!! (2/^A)'/^^Vi,e-'^^?^/^y,„(^i2, ^12). (33) 

The wave function is normalized according to 

J dr\<l>tiri2)f = 1. (34) 

A similar wave function (i>^^t{rzi) can be written for the relative motion between the quark q and the antiquark Q 
in the meson S(34). 

With the knowledge of the Gaussian wave functions and the effective interquark interaction, the meson-meson 
potential V(r) can be evaluated analytically. The meson-meson potential depends on {AxBxi\Vi\AqBq) which can be 
further decomposed as a sum of {A\B\i\v{rjk)\AoBQ), 

2 4 

{AxBx'\Vi\AoBo) =J2Y1 CjCk{AxBx'\v{rjk)\AoBo). (35) 

j = l k=3 

The matrix elements {AxBx'\v{rjk)\AoBo) can be evaluated in the configuration space where there are three indepen- 
dent coordinates because the coordinates of the four particles are connected by the center-of-mass coordinate ili234. 
We can fix the center-of-mass coordinate iii234 to be at the origin. 



-R1234 = (mi2-Ri2 + m34i?34)/(mi2 + 77234) = 0, 



(36) 
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and choose the independent coordinates to be r, ri2, and r34. Using Eqs. (2), (5), (6), (15), and (36), we can represent 
all other coordinates in terms of these three coordinates: 

Vjk = r + fA{jk)ri2 + fB{jk)r34, (37) 

where 

/.4(14)= /2, /b(14)= h, 

(38) 



/.4(14) = 


/2, 


/b(14) = 


/3, 


/a(13) = 


/2, 


/b(13) = 




/a (23) = 


-fl, 


/b(23) = 


-/4, 


/a(24) = 


-fl, 


/s(24) = 


/s- 



Using the above relations between rjk and {r, ri2, r34}, the matrix element (^A-BA'k('"jfc)|^o-Bo) can be written as 
{AxBx'\v{rjk)\AoBo) = j dr^ j dr34Pxo{ri2)px'o{r34)v{r + fA{jk)ri2 + fB{jk)r34), (39) 
where p^Q and pf,g are 

pUr 12) =rx{ri2)Mr 12), (40) 

PA'o(»-34) =<^A'(r34)<^0(r34). (41) 

Wc note that p^gC^is) i'' the density distribution of meson ^o(12) and Pxq{ti2) is the transition density for making the 
transition from the ground state y4.o(12) to the excited states ^4^(12). Similarly, p^o('''34) ^^"^ Pa'o(''"i2) '^^''^ respectively 
the density distribution and transition density of meson -B(12). By the method of Fourier transform [58, 68], the 
evaluation of the matrix element (39) can be greatly simplified. The result is 

{AxBy\v{rjk)\AoBo) = J ^^e^^'^ pio[fAUk)p] pUfB{jk)p] vip) (42) 



where 



pio^iP) = J dye^P Vpi^iy), (43) 



vip) = J drjke-'P-''^''virjk). (44) 

We first evaluate the matrix element {AoBo\v{rjk)\AoBo) which is also needed in subsequent calculations. From the 
wave functions we can determine 

p^fip) = e-P'/'^l^. (45) 
Therefore, the matrix element {AoBo\v{rjk)\AoBo) is related to v{p) by 

{AoBo\v{rjk)\AoBo) = J ^!|^e^P-^e-«>^t;(p) ^ voo{jk;a%,r), (46) 



where a^j^ is 



„2 _ [fAUk)r , [fBjjk)] 

- API + API ' 



and we have written (Ao-Bo|^^('''jfe)|^o^o) voo{jk;a'^f.,r) to exhibit fully its dependence on a^^, and r. 

The matrix element (y4.o-Bok('"jA;)|^o-So) can be expressed as the sum of contributions from the Yukawa interaction 
and the exponential interaction in Eq. (29), 

{AoBo\v{rjk)\AoBo) = (AoSo|«Y"'^(r,fe)|AoBo) + {AoBo\v'^P{rjk)\AoBo). (48) 
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The matrix element for the Yukawa interaction is found to be 

(AoBo|t;Y"^(r,fe)|AoBo) = ase'^'^'^^u^vi - U2V2), (49) 

ui = (1 + erf[(r - 2a%^l) / {2ajk)]) /2, (50) 

U2 = (1 - erf[(r + 2a%ii)l{2aiu)\) /2, (51) 

^^1 = , (52) 

r 

V2 = — . (53) 
r 

The matrix element for the exponential interaction is 

{AoBo\v''^'^{rjk)\AoBo) = — e"^^'* {mw^ + U2W2), (54) 

4/x 



where 



3 

and 



wi = (r - 20^^)— — , (55) 



W2 = {r + 2a%fi) — . (56) 

The direct potential vanishes when all the quark masses are equal, and it becomes more attractive as the mass of the 

heavy quark increases. 

We fix the vector r to lie along the z-axis, and quantize the azimuthal component of the angular momentum m 
to be projections along the z-axis. For the polarization potential V^oi('')) we consider intermediate excitations to the 
next excited states of ^(12) and -8(12) characterized by {A or \' = Im = \m}. It can be shown that the matrix 
element {A\Bxi\v{rjk)\AQBQ) is zero for {A = 0; A' = Ira] and {A = Im; A' — 0}. It gives a non- vanishing matrix 
element for {A = Im; A' = 1 — m} with {m = —1, 0, 1} when both mesons are excited. 

To evaluate the matrix element {AimB\-m\v{r jk)\AQBf)) , we write out the wave functions explicitly and obtain 

(AroB.oKr,.)|Ao5o) = / ?^^^e^^- e-^%^\\^v[p). (57) 

A comparison of the above equation with Eq. (46) yields 

{A,oB,oHrik)\A,B,) = -l^^^MJll^y^,^jk.,a%,r). (58) 
To evaluate the matrix element {AiiBi-i\v{rjk)\AoBo), we write out the wave functions explicitly and obtain 

{AnBi_i\v{rjk)\AoBo) = J ~(2wf~^ ^ 2 ^ ^ 

A comparison of the above equation with Eq. (46) gives 



With the analytical results of Eqs. (48), (49), and (54) for {AoBo\v{rjk)\AoBo) (or voo{jk;a'ji.,r)), the differentia- 
tion of i;oo(jfc; a|fe, r) with respect to a|j, and r can be readily carried out analytically to give the matrix elements 
{AimBi_rn\v{rjk)\AoBo) for {m = -1,0,1}. The energy denominator in Eq. (22), {e(Ai„)-|-e(i?i_„)-e(Ao)-e(i?o)}, 
can be obtained from experimental masses of D{1P), D*{1P), D{1S), and D*{1S). The masses of the IP states of 
the open bottom mesons are not known experimentally. We infer the mass differences from theoretical masses of open 
bottom mesons calculated in [42] where we get M{B(IP))- M{Bo) = 0.5348 GeV, and M(B*(1P))-M(B^) = 0.5168 
GeV. These quantities allows us to calculate the polarization potential using Eq. (22). 
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FIG. 3: The total potential V(r) and its various components for the DD* system. The direct potential Vdir is the sum of 
^^ukawa ^^exponential ^ V{r) is the sum of the direct potential and the polarization potential. The insert with an 
expanded ordinate scale gives a clearer picture of the total potential and its various components at large distances. 



V. HEAVY QUARK MOLECULAR STATES 

The meson wave function in Eq. (33) is characterized by a momentum width parameter /? that is related to the 
root-mean-square quark- antiquark separation frms of the meson by 

Previously, root-mean-square quark-antiquark separation of mesons have been obtained from calculations of meson 
wave functions and meson masses (Table IV of [42]). From this table, we use the root-mean-square quark-antiquark 
separations of rrms(-D) = 0.585 fm, r^ras{D*) = 0.626 fm, rrms(-B) = 0.574 fm, and rrms{D*) = 0.583 fm to calculate 
the /? values for various meson wave functions. In the present numerical work, we shall include only the screened 
color-Coulomb interaction (Yukawa interaction) and the screened linear interaction (exponential interaction) so as to 
study the gross features of the meson-meson potential and molecular states. 

Using the effective interaction outlined above, we obtained meson-meson potentials for various combinations of 
D,D*,B,B* with D, D* , B, B* . The total potential V{r) and its various components for DD* are shown in Fig. 3. 
The insert in Fig. 3 gives a more expanded description of various components at large distances. The total direct 
potential is the sum of the contributions from the Yukawa and the exponential interaction, and is attractive. The 
magnitude of the polarization potential is smaller than the magnitude of the direct potential at short distances but 
the role is reversed at intermediate separations. The two different contributions are about equal at large separations 
beyond 1.2 fm. As a result, the total potential, which is the sum of the direct potential and the polarization potential, 
has an extended attractive region extending to 1 fm and beyond. 

With such a potential and the experimental masses of D and D*, we solve for the lowest eigenstate of the meson- 
meson system. We find a bound state at a binding energy of 7.51 MeV. We plot the wave function for this molecular 
state in Fig. 4. The wave function ipir) is quite extended and has a significant amplitude beyond 1 fm. This state 
is weakly bound, and has a large root-mean-square D-D* separation of 1.37 fm. Its weak binding and its extended 
spatial separation between D and D* provide a good characterization of its molecular structure. 
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r (fm) 

FIG. 4: The wave function il){r) of the DD* system as a function of the separation between D and D* . 

For a flavor-independent effective interaction, the theoretical meson-meson potentials and the binding energies are 
expected to change very insignificantly for different flavors of the hght quarks. On the other hand, the mass of the 
system is related to the binding energy B by 

M = M(A,) + M{B,,) + e = M{A,) + M{B,,) - B, (62) 

and the masses of and B^i depend on light quark flavors. Therefore, the calculated mass of the {Qq)-{qQ) molecular 
system depends on light quark flavors. We list the molecular states for different combinations of the {cq) and {qc) 
mesons in Table I. The properties of the states in the DD, DD*, and D*D* systems are quite similar. 

Table I. The calculated masses, binding energies, and the root-mean-square {cq)-{qc) separation rrms of molecular 
states formed by (cq) and (qc) mc^sons. 



{cq)-{qc) 


State 


Binding Energy (MeV) 


Mass (MeV) 


rrms (fm) 


D'^D'^ 


IS 


3.10 


3725.90 


1.73 


D^D-, D+D° 


IS 


3.10 


3730.70 


1.73 


D+D- 


IS 


3.10 


3735.50 


1.73 


D°D*° 


IS 


7.53 


3863.67 


1.37 


D^D*- 


IS 


7.53 


3866.97 


1.37 


D+D*° 


IS 


7.53 


3868.47 


1.37 


D+D*- 


IS 


7.53 


3871.77 


1.37 




IS 


15.31 


3998.09 


1.06 




IS 


15.31 


4001.39 


1.06 


D*+D*- 


IS 


15.31 


4004.69 


1.06 



To produce the unknown state X(3872) in the process B^ K^X{2>&T2) observed by the Belle Collaboration, the 
state X(3872) most likely originates from the production of a cc pair with / = during the weak decay of the bottom 
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quark. On the other hand, the calculated masses of D^D*^ and D~^D*~ molecular systems are respectively 3863.67 
MeV and 3871.77 MeV (see Table I). By the symmetry of charge conjugation, the masses of D*^D^ and D-D*+ 
systems arc also 3863.67 McV and 3871.77 McV respectively. The most suitable candidate for the 3872 McV state 
appears to be the molecular state D^D*^ and its charge conjugate D~D*~^. One expects that the molecular states 
D^D*^, D~^D*~ , and D~D*~^ will be mixed to form components of 7 = and 1 = 1 states. As they lie in the vicinity 
of the 3872 state, it is tempting to identify the 7 = of this set of multiplets as the 3872 state observed by the Belle 
Collaboration. It has been conjectured however that the 3872 state may break isospin symmetry maximally [4, 6] and 
more experimental information are needed to understand better the isospin properties of the 3872 state. 
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B+B* Potential 
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.0 0.2 0.4 0.6 0.8 1.0 1.2 

r (fin) 



FIG. 5: The total potential V{r) between B and B* and its various components, 
various components at large separations in greater detail. 



The insert gives the total potential and its 



We show in Fig. 5 the potential for the BB* system. The potential has features similar to those of the DD* 
system except that the BB* potential is deeper than the DD* potential. This arises because the quark mass ratio 
fnb/mu,d is greater than mc/mu,d- For the BB* system, the magnitude of the direct potential is consistently much 
larger than the magnitude of the polarization potential. The total potential is attractive at large separations. We 
search for eigenstates for the BB* system. The large mass of the bottom quark and the large depth of the BB* 
potential leads to two eigenstates of the BB* system. There is a IS state with a binding energy of 151.75 McV and 
a root-mean-square B-B* separation of 0.314 fm. In addition, there is a 2S state with a binding energy of 0.88 MeV 
and a root-mean-square B-B* separation of 2.23 fm. 

The wave functions for these two states arc shown in Fig. 6. The ipis{r) wave function is confined at short distances, 
but the ip2S {r) wave function is spatially quite extended. We list in Table II the calculated binding energies, masses, 
and the root-mean-square radii of states formed by (bq) and {qb). As the masses of the B^ and B"^ are nearly the 
same, we do not need to distinguish B (and similarly B*) mesons of different light quark flavors. 

Table II. The calculated masses, binding energies, and the root-mean-square {bq)-{qb) separation rrms of molecular 
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states formed by (bq) and {qb) mesons. 



m-{qb) 


State 


Binding Energy (MeV) 


Mass (McV) 


Trrns (fm) 


BB 


IS 


148.24 


10410.56 


0.317 


BB* 


IS 


151.75 


10452.64 


0.314 


B*B* 


IS 


155.40 


10494.60 


0.311 


BB 


2S 


0.80 


10558.00 


2.25 


BB* 


2S 


0.88 


10603.52 


2.23 


B*B* 


2S 


0.95 


10649.05 


2.21 



20.0 - 



B-B* wave function 




■5. 



to 



J_ 



J I L_ 



J_ 



1.0 



3.0 



2.0 

r (fm) 

FIG. 6: The wave function tp{r) of the BB* system as a function of the separation between B and B* 



The potentials and wave functions for {cq)-{qb) systems involving both charm quark and bottom quarks are similar 
to those of BB* and DD* and have similar IS* eigenstates with a binding energy about 29-43 MeV and a root-mean- 
square {cq)-{qb) separation about 0.6-0.8 fm. We list in Table III the eigenstates for various {cq)-{qb) systems. 



Table III. The calculated masses, binding energies, and the root-mean-square {cq)-{qb) separation rr, 
states formed by (eg) and (^6) mesons. 



of molecular 



{cq)-{qb) 


State 


Binding Energy (MeV) 


Mass (MeV) 


Trms (fm) 


D^'B 


IS 


27.83 


7115.01 


0.756 


D+B 


IS 


27.83 


7119.81 


0.756 


D'^B* 


IS 


29.36 


7158.96 


0.739 


D+B* 


IS 


29.36 


7163.76 


0.739 


D*^B 


IS 


40.48 


7244.07 


0.638 


D*+B 


IS 


40.48 


7247.38 


0.638 


D*^B* 


IS 


42.63 


7287.34 


0.624 


D*+B* 


IS 


42.63 


7290.64 


0.624 
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VI. POTENTIAL FOR THE {Qq)-{Qq) SYSTEM 

From our present model of effective charges, it is easy to see how the direct potential for the {Qq)-{qQ) system 
depends on the ratio of the heavy quark mass mg to the light quark mass m,. If the mass of Q is the same as the mass 
of q, then the attraction between effective charges of opposite signs will be compensated by the repulsive interaction 
between effective charges of the same sign, when we bring the two mesons together. The net direct potential is zero 
when mg = rUq. The polarization potential is always attractive, but the magnitude of the polarization potential in 
most cases is not attractive enough to lead to a bound meson-meson system. 



0.8 



0.6 



^ 0.4 
S 0.2 



0.0 



-0.2 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

r (fm) 

FIG. 7: The total potential V{r) between D and D is shown as the solid curve. Shown also are various components of the 
potential. The insert gives the potentials at large separations in greater detail. 

If mg 3> TOg, the light quark and antiquark in {Qq)-{qQ) will move around a pair of nearly static heavy quark and 
heavy antiquark. As one brings a {Qq) in the vicinity of a {qQ) meson, the light quark and the light antiquark will 
attract each other because of their opposite effective charges and the meson-meson interaction will be attractive. The 
greater the mass ratio mg/mq, the greater their attraction will be. Furthermore, the attractive polarization potential 
extends to large separations and facilitates the formation of molecular states. 

From the above argument, one notes that the direct potential will be repulsive for the {Qq)-{Qq) system which 
has two heavy quarks and two light antiquarks. In this case, as the mesons are brought close together, the two 
antiquarks moving around the nearly static heavy quarks will interact. As the effective charges of the antiquarks have 
the same sign, they tend to repel each other and the resultant direct interaction is repulsive. We show in Fig. 7 the 
total potential and its various components for the DD system which has a repulsive direct potential. Although the 
polarization potential is attractive, the total potential is only very weakly attractive between 0.5 fm and 1 fm with 
a maximum depth of about 6 MeV. Such a shallow potential is not attractive enough to hold a molecular state. It 
should be pointed out, however, that symmetrization or anti-symmetrization of the spatial wave functions of the light 
antiquarks orbiting around the heavy quarks will give rise to additional effects [36, 37]. These effects must be studied 
in more detail for the QQqq system, as in the case of electrons in a diatomic molecule. In these future studies, the 
direct potential and the polarization potential will provide some of the information needed for such an investigation. 
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VII. CONCLUSION AND DISCUSSIONS 

We have studied four-quark {Qq)-{qQ) systems involving charm and bottom quarks in a quark-based model. The 
system is described by a non-relativistic four-body Hamiltonian with a pairwise interaction. 

We note that the molecular state of our interest is weakly-bound and the average separation R between the heavy 
quark mesons is considerably greater than the average radius a of the mesons. An R » a, the dominant process 
for the interaction of the mesons is the exchange of two gluons, which leads to the color van der Waals interaction 
[53-55]. We can equivalently represent the color van der Waals interaction in terms of effective charges for quarks 
and antiquarks in a QED-typc interaction. These effective charges give the proper description of bound states for 
a quark interacting with an antiquark in an isolated mesons. They also gives rise to a color-clcctric dipole-dipole 
interaction between mesons when a meson is brought in close proximity to another meson [53 55] . The color van der 
Waals interaction should be screened by introducing a screening mass /i due to the breaking of the gluonic string 
at large distance. Incorporating these concepts in a phenomenological interquark potential, we find a weakly-bound 
£)+£)*- state near the threshold that may be qualitatively identified as the 3872 state observed recently by the Belle 
Collaboration. 

In the present model, the meson-meson potential for the {Qq)-{qQ) system can be decomposed into a direct potential 
and a polarization potential. The direct potential depends on the ratio of the quark masses mg/mq. If the mass 
ratio is 1, then the direct potential is zero. The magnitude of the attractive direct potential increases as the mass 
ratio increases. The polarization potential is always attractive. However, only when the quark mass ratio of both 
mesons arc sufficiently high can the total potential be attractive enough to lead to bound states. There are no 
bound molecular states for {sq)-{qs) and {cq)-{qs) systems if we use meson-meson potentials obtained in the present 
model. We find bound meson-meson states in open heavy-quark systems involving combinations of {D, D* , B, B*} 
with {D, D*,B, B*}. For open bottom meson pairs, we found weakly-bound 2S states near the threshold, in addition 
to deeper IS states with binding energies of about 150 MeV. 

We have included only screened color-Coulomb and screened linear confining interactions in the present work. It 
will be of interest to study in future work the influence of the spin-spin and other components of the interaction on 
the molecular structure. In meson systems having spin zero and one, the spin-spin interaction does not contribute to 
the direct potential when the spin of one of the mesons is zero. The effect of spin-spin interaction is expected to be 
small when one of the mesons is D, B, D, or B. 

The present model differs from those of the pion-exchange models and other previous multi-quark models [3- 
9, 17-41, 45-49]. The results of molecular states obtained from different models can naturally be different. The 
present model predicts many molecular states in the combinations of {D, D* , B, B*} with {D, D* , B, B*}. While it 
is encouraging that the present phenomenological model gives a molecular D~^D*~ state at about the right energy 
as the 3872 state, further careful comparisons of different predictions with experiment will be needed to determine 
whether the 3782 state is a molecular state. Close and Page [4] and Voloshin [6] suggest that if the 3872 state is 
a molecular state, its D and D* components should decay with a width equal to those of the isolated mesons, and 
the 3872 state should also be seen in the invariant mass of of the combination of the decay products of D and D* . 
Furthermore, while the present model gives DD and BB molecular states, the pion-exchange model [3, 23] does not 
predict DD and BB molecular states. A search for DD and BB moleciflar states will be of interest to distinguish 
the different mechanisms that are present in generating the molecular states. 
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